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ABSTRACT 

We analyze the distinguishability of populations of coalescing binary neutron stars, 
neutron-star black-hole binaries, and binary black holes, whose gravitational-wave sig¬ 
natures are expected to be observed by the advanced network of ground-based inter¬ 
ferometers LIGO and Virgo. We consider population-synthesis predictions for plau¬ 
sible merging binary distributions in mass space, along with measurement accuracy 
estimates from the main gravitational-wave parameter-estimation pipeline. We find 
that for our model compact-object binary mass distribution, we can always distin¬ 
guish binary neutron stars and black-hole-neutron-star binaries, but not necessarily 
black-hole-neutron-star binaries and binary black holes; however, with a few tens of 
detections, we can accurately identify the three subpopulations and measure their 
respective rates. 
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1 INTRODUCTION 


Merging binary neutron stars (NS-NS), binary black holes 
(BH-BH), and mixed black-hole - neutron-star binaries 
(BH-NS) are among the most likely sources for the ad¬ 
vanced gravitational-wave detectors LIGO (Harry, G. M. & 


the LIGO Scientihc Collaboration 20101 and Virgo (Virgo 
Collaboration|2d09 ). The rates could range from a detection 
every few years to a few hundred detections per year for each 


source type (Abadie et al. 20101 once detectors reach full 


sensitivity toward the end of the decade, although the hrst 
detections are possible as early as the end of 2015 (LIGO 


Scientihc Collaboration et al. 20131. One key question is 


whether it will be possible to determine the type of the bi¬ 
nary — NS-NS, BH-NS, or BH-BH — from its gravitational- 
wave signature. 

The surest way of distinguishing NS and BH binaries is 
via tidal dissipation in neutron stars leading to additional 
loss of orbital energy on top of that lost through gravita¬ 
tional waves (GWs), and their possible eventual tidal dis¬ 
ruption. However, the dominant GW signature through most 
of the LIGO-Virgo frequency band will be well-described by 
point-particle waveforms, with tidal effects becoming impor- 
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tant only at higher frequencies where detectors are less sen¬ 
sitive, making it difficult to distinguish binary types in this 


way except, perhaps, for the loudest sources (e.g.. Read et al. 
|2013[ and references therein). Electromagnetic signatures in 
the form of short gamma-ray bursts, afterglows, and kilono- 
vae could also help distinguish binaries with matter that 
can undergo tidal disruption (NS-NS or BH-NS) from those 
without (BH-BH); however, electromagnetic signatures may 
prove difficult to observe even for most NS-NS and BH-NS 
binary mergers (e.g., [Kelley et al.|2013p , nor could they dis¬ 
tinguish between those binary types given the current lim¬ 
ited understanding of emission processes. Other approaches 
to distinguishing subpopulations of different binary types in 
the overall set of GW observations would rely on classifica¬ 
tion based on clustering in the parameter space of all avail¬ 
able information about the binaries, including spin, which 
appears to be observationally higher for BHs than NSs in 
compact binaries (e.g., Mandel fc O’Shaughnessy|2010 and 
references therein). However, in the absence of confident pre¬ 
dictions for NS and BH spin, here we will focus on classifi¬ 
cation based exclusively on component mass. 


Mass-based binary classification depends on (i) the as¬ 
sumptions about the mass population of NS and BH bina¬ 
ries, (ii) the accuracy with which component masses can be 
determined from the GW signature of an individual merger 
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observation, and (iii) inference about the source population 
based on a statistical analysis of multiple detections. In this 
Letter, we address all three of these issues in turn. We argue 
that (i) astrophysical models predict that NS-NS, BH-NS, 
and BH-BH binary subpopulations may be strongly clus¬ 
tered, aiding source identification; (ii) a careful analysis of 
parameter measurement accuracy is required to take into 
account possible degeneracies; (iii) even if the binary type 
cannot be confidently determined from an individual mea¬ 
surement, a set of observations can be used to constrain 
mass distributions and measure the rates of different merger 
types. We find that for our model of the compact-object bi¬ 
nary mass distribution, based on population-synthesis sim¬ 
ulations (Dominik et al. 20141, we can always distinguish 
NS-NS and BH-NS binaries, but not necessarily BH-NS and 
BH-BH binaries. However, with a few tens of detections, 
we can identify subpopulations and measure their respec¬ 
tive rates to an accuracy dominated by the Poisson count¬ 
ing statistics of each binary type rather than contamination 
among subpopulations. 


2 MASS DISTRIBUTIONS 


If the NS and BH mass distributions are directly adjacent, 
then it will be impossible to conhdently determine the type 
of an individual source from finite-accuracy mass measure¬ 
ments alone. However, observations of NS binaries, particu¬ 
larly via radio pulsars, and BH X-ray binaries indicate that 
there may be a mass gap between the two populations: the 


highest measured NS masses just exceed 2Mq (Lattimer 


20121 while BH masses may only start at ~ 4 - 5. 5Mf? 


depending on the assumed shape of the distribution (Ozel 
et al.|20T0 Farr et al.|2011 1; but see Kreidberg et al. (20121 


for possible selection biases that could push the BH masses 
lower. This mass gap could be indicative of the supernova 
explosion mechanism (Belczynski et al. 

Previous analyses of the distinguishability of source 
populations based on GW observations typically relied on 
a mass gap individually applied to both components (e.g., 
Hannam et al.|[2013l. However, in practice, even if NSs oc¬ 


cupy the mass space through 2 solar masses and BH masses 
start at 5 solar masses, it is not necessarily the case that 
merging compact binaries can contain arbitrary combina¬ 
tions of NS and BH masses within those respective ranges. 
The compact-object mass distribution in binaries depends 
on the outcome of binary evolution. To carry out a self- 
consistent analysis of the distinguishability of source popu¬ 
lations, we consider sources that arise from an astrophysi- 
cally modeled population of binaries and compare their mass 
measurements against the joint boundaries on component 
masses as given by the same population. 

For this study, the binary population was generated 


with the StarTrack population synthesis code (Belczynski 
|et al.|200^ under the “Stand ard” model B of|Dominik et al. 
(20121. This model builds on Belczynski et al. (20081 by in¬ 


cluding the rapid supernova engine of Fryer et al. (20121, 


which allows for the formation of NSs with large natal kicks 
from single stars with zero-age main sequence (ZAMS) mass 
below ~ 2 IM 0 , and nearly complete fallback onto BHs with 
correspondingly low natal kicks from more massive progen¬ 
itors, yielding minimum BH mass above ~ 5Mq; realis¬ 
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Figure 1. The mass distribution in Me — rj space of compact- 
object binaries detectable by advanced gravitational-wave obser¬ 
vatories a s simulated by ||Dominik et al. |2014| l, using the standard 
model of | |Dominik et al.|2012[ |. NS-NS, BH-NS, and BH-BH pop¬ 
ulations occupy the top left (red), bottom middle (green) and 
top right (blue) regions of the plot, respectively. Each dot cor¬ 
responds to a potentially detectable source, but the number of 
binaries plotted is increased far beyond the anticipated detection 
rate | |Abadie et al.|[2010^ to better indicate the population dis¬ 
tribution. Chirp mass boundaries for the three distributions are 
shown with dashed red, dash-dotted green, and dotted blue ver¬ 
tical lines, respectively. BH-BHs with chirp mass above 20Mq, 
which cannot be confused for another source type, are omitted 
from the figure. 


tic prescriptions for common-envelope binding energies and 
mass loss through stellar winds; potential binary disruptions 
in supernovae and stellar mergers within the common en¬ 
velope; and a prescription for electron capture supernovae 
which lead ~ 7-11Mq ZAMS mass stars in binaries to form 


NSs with low natal kicks and masses near I. 25 M 0 ; see (Do¬ 


minik et al. 20121 for details. We note that many of the input 


parameters and assumptions that go into population synthe¬ 
sis, such as those related to mass loss and mass transfer, BH 
supernovae fallback and kicks, and common envelope evo¬ 
lution, are poorly constrained. Therefore, this is only one 


of many possible model variations; see the series (Dominik 
et aI.|20T2 2013 20141 for discussion, and (Postnov & Yun 


gelson||2014 1 for a recent review of decades of work in this 


held. We do not insist on this model being accurate, but only 
use it for illustration purposes to show how general features 
of the model, and particularly clustering of the subpopula¬ 
tions in mass space, can be used to distinguish the compact 
binary type. 

Figure shows the distribution of masses for the three 
binary types among the binaries detectable by the advanced- 
detector network as estimated by ( Dominik et al.|[2014 |. A 
wide range of metallicities is included in these calculations 
to simulate the redshift-dependent metallicity distribution 
in the Universe; in fact, as (Dominik et al. 20141 showed, 
the local BH-BH mergers are dominated by low-metallicity 
massive binaries that formed in the early Universe, with long 
time delays between star formation and merger. 

We re-parametrize the component masses mi, m 2 via 
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the chirp mass Me = (rn\ + , which is 

generally well-measured with gravitational waves, and the 
symmetric mass ratio rj = mim2(mi -|- m 2 )~^. The main 
features of the model become apparent from the figure: the 
three source types are in general clearly separated in the 
chirp mass — mass ratio space. 

NS-NS systems are very strongly clustered near equal 
masses (77 > 0.2477, or m 2 /mi > 0.82). This clustering of 
the NS-NS mass ratio reflects their nearly equal- (and typ¬ 
ically low-) mass progenitor binaries, preferred since they 
allow the primary to form through an electron-capture su¬ 
pernova which avoids disrupting the binary; mass transfer 
after the hrst supernova leads to a growth in the mass of 
the primary by no more than ~ 7%, keeping the mass ratio 
nearly equal. BH-BH systems also cluster near equal masses 
at solar metallicity, in part because a large difference in ini¬ 
tial masses causes unstable mass transfer and a common- 
envelope merger during the first mass transfer phase, and 
in part because of the very efficient mass loss through stel¬ 
lar winds which tends to cluster BH masses despite a very 
broad range of ZAMS masses. Mass loss rates are reduced 
for lower metallicities; therefore, more unequal mass ratios 
are possible for lower metallicities, with mass ratios down 
to ?7 ~ 0.1 for BH-BHs that formed at 0.01^©. Meanwhile, 
BH-NS systems have more extreme mass ratios, ranging be¬ 
tween rj ~ 0.03 and rj ~ 0.2, consistent with the supernova 
fallback prescription that naturally generates the mass gap. 

The chirp mass is lowest for NS-NS systems given 
the constraints on the maximum component masses. In 
fact, there’s a gap between the highest NS-NS chirp mass 
(Me = 1.67M0, for a 1.92M0 -|- I.9IM0 binary formed at 
O.IZ0) and the lowest BH-NS chirp mass (Me = 2.IIM0, 
for a 6.IOM0 -I- I.IOM0 binary formed at O.75Z0). Mean¬ 
while, BH-BH systems do overlap with BH-NS binaries in 
chirp mass, as the lowest-mass equal-mass BH-BH systems 
(5.55M0 -|- 5.53M0 forming at O.O25Z0) can have a lower 
chirp mass (Me = 4.82M0) than the highest-chirp mass, 
extreme-mass-ratio BH-NS system (Me = 6.42M0 for a 
48.OM0 -I- I.7OM0 binary that formed at 0.005Zq). 


3 MEASUREMENT ACCURACY 


We estimate the accuracy of mass measurements on individ¬ 
ual events through the LALInference parameter-estimation 
toolkit, designed for efficient stochastic exploration of the 
parameter space of GW signals ( Veitch et al.|[2015 l. We es¬ 
chew the large-scale mock data challenges attempted else¬ 
where (e.g., Singer et al.|2014 Berry et al. 20141 and instead 
consider several events placed at the boundaries of the re¬ 
gions described in the previous section. 

We analyze mock injections of the highest-chirp-mass 
NS-NS system labeled as detectable in our population- 
synthesis simulation, with (1.92M0, I.9IM0) compo¬ 
nents; the lowest-chirp-mass BH-NS system with (6.IOM0, 
I.IOM0) components; the highest-mass-ratio BH-NS system 
with (5.6OM0, 2.IOM0) components; the highest-chirp-mass 
BH-NS system with (48.OM0, I.7OM0) components; and the 
lowest-chrip-mass BH-BH system with (5.55M0, 5.53Mq) 
components. Given the uncertainty about NS and BH spins, 
we consider four variations of each injection: non-spinning 
components, components with aligned spins, and spinning 



Figure 2. Posterior PDF samples for the 20 injections de¬ 
scribed in the text: NS-NS (red), BH-BH (blue), and three BH-NS 
(green), each with four spin variants. Point opacity corresponds 
to the relative value of the posterior probability density function. 
Chirp mass measurements are typically very accurate (with the 
exception of the very massive, extreme-mass-ratio BH-NS binary 
on the bottom right of the plot), while mass ratios are partly de¬ 
generate with spins, and their measurement accuracy depends on 
the exact injected spin configuration. 


components with two randomly chosen arbitrary spin direc¬ 
tions. In all cases when components are spinning, the BH 
dimensionless spin parameter is y = 0.8; the NS is not spin¬ 
ning in BH-NS injections and only the higher-mass NS is 
spinning with a dimensionless spin parameter of x = 0-05 
in spinning NS-NS injections (representative of the fastest 
known NS spin in a confirmed merging Galactic NS-NS bi¬ 
nary, J0737-3039A). For NS-NS and BH-BH systems, as well 
as the highest-mass-ratio BH-NS system, we use for all injec¬ 
tions and templates the IMRPhenomP precessing waveform 
family ( Hannam et al.|2014 1, which models all phases of the 
waveform: inspiral, merger, and ringdown. However, since 
this waveform model is restricted to mass ratios less extreme 
than 1 : 4, i.e., rj 0.16, BH-NS injections (other than the 
high-mass-ratio injection) are made and analyzed with Spin- 
TaylorT4 waveforms that include only the post-Newtonian 
inspiral phase ( Buonanno et al?]|2003 ), not the merger and 
ringdown (see Littenberg et al.|2015 for an analysis of BH- 
BH signals with SpinTaylor waveforms, which shows some¬ 
what greater typical mass measurement uncertainty). 

In all analyses, we allow for arbitrary, precessing spins in 
the templates. We use flat priors on the component masses, 
flat priors on dimensionless spin magnitudes x [0i 1] 
(x [0, 0.9] for IMRPhenomP analyses because of limita¬ 
tions in that waveform’s region of validity) and isotropic 
priors on spin directions for the Bayesian analysis. All anal¬ 
yses are carried out from a starting frequency of 40 Hz with 
mock signal-only data from a LIGO-Virgo detector net¬ 
work operating at advanced LIGO sensitivity in the zero- 


detuned, high-power conhguration ( LIGO Scientific Gollab- 
|oration||2010[ ). The injected distance was chosen to yield a 
network signal-to-noise ratio of 12, roughly at the lower limit 


of detectability (e.g., LIGO Scientific Collaboration et al. 


|2013[ [Berry et al.||2014[ ), thereby yielding conservative pre- 

dictions on parameter-estimation accuracy. 
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Figure shows the posterior probability density func¬ 
tions (PDF) for the 20 analyses described above (four injec¬ 
tions - non-spinning, spin-aligned, and two processing - for 
each of 5 events). Several thousand samples for each event 
are drawn from the posterior by LALInf erence; the opacity 
of plotted points rises with the value of the posterior PDF. 
The injected masses are denoted with black crosses. 

The chirp masses are typically very well measured (pos¬ 
teriors look like narrow nearly-vertical lines), except for 
the extreme-mass-ratio, high-mass BH-NS injection that ap¬ 
pears on the bottom right of the plot. This is consistent 
with other analyses (e.g., Aasi et n(]|2013 Rodriguez et al. 


|2014[ [Hannam et al.||20l3 l, and is readily understood from 

the governing role played by the chirp mass in the phase 


evolution of the binary during the inspiral (e.g., Buonanno 
|et al.|[^09[ |Ohme et al.||2013[ ). The 90% credible regions in 
chirp mass spanning all four spin variations for the given 
component mass combinations are < O.OIMq for the NS-NS 
injections and < O.IMq for the BH-BH injections. For BH- 
NS injections, the 90% credible regions in chirp mass span 
< O.OSMq, except for the very poorly measured extreme- 
mass-ratio, high-mass injections discussed below, for which 
the combined 90% credible region in chirp mass spans ~ 
2Mq. Hence, measurement of the chirp mass alone allows us 
to distinguish NS-NS and BH-NS binaries, which are sepa¬ 
rated by a chirp-mass gap of O.45M0 - much larger than the 
chirp-mass measurement uncertainty - if we trust the model 
described in the previous section. Conversely, a single source 
falling into the chirp-mass gap between NS-NS and BH-NS 
binaries would allow us to disprove this model. Of course, 
we do not necessarily trust the population-synthesis predic¬ 
tions for the exact boundaries of the NS-NS, BH-NS, and 
BH-BH clusters in the mass plane; this will be addressed in 
the following section. 

On the other hand, the mass ratio r] is more poorly 
measured, and partly degenerate with spin (e.g.,|Ohme et ah 


|2013[ ). The difficulty of measuring the mass ratio is apparent 
from the very elongated posteriors in Fig. and from the 
typical “banana-shaped” plots of posteriors in component- 
mass space shown in Fig. The latter plot moreover illus¬ 
trates the dependence of measurement accuracy on the con¬ 
figuration of spins in the injected system. Systems with pre¬ 
cession typically have the best-measured masses at a fixed 
signal-to-noise ratio, as evidence of precession allows for 
more accurate spin measurements, partly breaking the mass- 
ratio-spin degeneracy ([van der Sluys et al.|2008||Chatziioan-| 


nou et al. 2014 20151. The limited accuracy of measuring the 


mass ratio rj means that posteriors for BH-BH and BH-NS 
events could overlap, making it difficult to distinguish be¬ 
tween the two. However, our present attempts to estimate 
the extent of the posterior for massive, extreme-mass-ratio 
BH-NS binaries that could overlap with the BH-BH distri¬ 
bution are limited by the lack of waveforms that include 
both precession and mergers and ringdowns for such sys¬ 
tems. By injecting and analyzing the BH-NS system on the 
bottom right of Fig. which has a total mass of ~ 50Mq, 
with inspiral-only waveforms, we are losing significant infor¬ 
mation and may be artificially extending the measurement 
uncertainty, creating, in our plot color scheme, “the green- 
eyed monster which doth mock”. Moreover, if we are suffi¬ 
ciently fortunate to have significant precession in the bina¬ 
ries Chatziioannou et al. (see Fig. [^and| |2015[), the posteriors 



Figure 3. 90% credible regions for the four BH-BH injections in 
component mass space (mi defined to be the larger component 
mass, mi > m 2 ). All injections have 5.55Mq, 5.53Mq compo¬ 
nents, but differ in spin magnitudes (0 for NoS, 0.8 for each of 
the dimensionless spins for the other injections) and directions 
(aligned for AS, random and processing for the two PS injections). 
Mass measurements are more accurate when precession breaks 
degeneracies between intrinsic parameters. The second process¬ 
ing injection may be atypically favorable for mass measurement 
(cf. [Littenberg et ar]|2015[ l and is shown as an example of what 
could be achieved for a nearly edge-on processing binary. 


will be narrowed and more accurate parameter estimation 
will be possible. Furthermore, we considered parameter es¬ 
timation at a signal-to-noise ratio near the threshold for de¬ 


tectability (e.g., LIGO Scientific Collaboration et al. 2013 


Berry et al.|2014 l; louder signals will lead to more accurate 

inference. 

In any case, if the population boundaries are known in 
advance (e.g., if we have sufficient trust in population syn¬ 
thesis models models of Sec. we anticipate being able to 
readily distinguish S> 99% of the BH-NS and BH-BH popu¬ 
lations on the basis of chirp mass measurements alone for the 
model we have considered: the 99th percentile of the BH-NS 
chirp-mass distribution, 4.76M0, is readily distinguishable 
from the 1st percentile of the BH-BH chirp-mass distribu¬ 
tion, 5.67M0, given chirp mass measurement uncertainties 
< O.1M0. 


4 POPULATION STATISTICS 

We now focus on the question of distinguishing types of com¬ 
pact binary sources within the observed population. Given 
the uncertainty in the claimed range of a putative mass gap 
between NS and BH masses, and the fact that the mass 
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combinations of merging binaries may be more restrictive 
than suggested by the mass gap, as discussed in Sec. 
we address the question of classification using only the ob¬ 
served population itself, rather than relying on any a priori 
assumptions about subpopulation boundaries or mass gaps. 
We provide approximate scaling arguments and typical use¬ 
ful constraints. The exact accuracy of classification depends 
on the details of the true distribution; the plausible two- 
dimensional mass distribution described in Sec. [2] is used as 
an illustration. 


Measuring a distribution boundary from the observed 
population. A boundary of a distribution, if abrupt, is gen¬ 
erally easier to measure than other distribution parameters 
(cf. Mandel et al. J2^. For example, with perfect accuracy 
on individual measurements, the accuracy of the estimate 
of the mean of a sampled distribution scales inversely with 
the square root of the number of observations (samples) - 
but the location of a sharp boundary in the distribution can 
be estimated with an accuracy that scales inversely with 
the number of observations. If observations have an indi¬ 
vidual measurement uncertainty a, and N observations are 
spread over a parameter range A, an abrupt boundary can 

be estimated with accuracy ~ \/7^ + course, the 

boundary is harder to measure if the edge is gradual rather 
than abrupt. However, our specific two-dimensional mass 
distributions do, in fact, have fairly sharp one-dimensional 
boundaries following suitable coordinate transformations; 
e.g., in our model, the mass of the lower-mass BH-BH com¬ 
panion has a boundary of m2 > 5.5Mq. The smaller com¬ 
panion mass has a one-sigma uncertainty of ~ I.2M0 for 
comparable-mass binaries (taking the worst-case credible re¬ 
gion of Fig. which shows 90% credible intervals), and the 
BH-BH m2 population width is A ~ IOM0, so ~ 10 ob¬ 
servations conhdently known to be BH-BHs are needed to 
constrain the lower boundary on BH mass in the BH-BH 
distribution to an accuracy of < IM0. 


Establishing the existenee of a break in the distribution. 
Observing the existence of a break between subpopulations 
need not rely on very precise individual measurements. In 
fact, even perfectly accurate measurements cannot individ¬ 
ually confirm the presence of break in the distribution (al¬ 
though a single perfect measurement that falls into an ex¬ 
pected break would disprove its existence). Rather, the ex¬ 
istence of a break is indicated by a drop in the density of 
observations in the parameter range of the break. Critically, 
this is observable even when the parameter measurement 
accuracy is comparable to or even exceeds the width of the 
break. Denoting the width of the break by Abreak (as op¬ 
posed to the net distribution width A), the measurement 
error in the parameter direction across the break by a, and 
the total number of observations by N, the expected num¬ 
ber of observations in the break is ~ erf ^ j. By 
comparison, if there were no break in the underlying dis¬ 
tribution, A'^Abuak observations would be expected in that 
region. Therefore, the dearth in the number of observations 
in the break will be significant when the difference between 
the actual number of observations and N jg large rel¬ 
ative to the Poisson fluctuations in the latter, N ^ 

[This does not indicate an absolute break, which cannot be 
established with a hnite number of observations, but a sig- 


nihcant local drop in the probability density function of the 
inferred underlying distribution.] Again, using the break be¬ 
tween BH-NS and BH-BH distributions as an example, with 
Abreak ~ 3.3Mq in m2, approximately 60 observations are 
required for a confident detection of a break at a three-sigma 
signihcance. In a contemporaneous study, [Littenberg et alT] 
(20151 find that hundreds of observations may be necessary 
to measure a break or gap in the mass distribution when 
considering an ad hoc population that is flat in component 
masses; this is consistent with the application of our analysis 
to their distribution. 

Classifying observations into sub-populations. As men¬ 
tioned earlier, we do not wish to rely on a priori divisions of 
the compact-binary mass parameter space into regions, but 
instead to measure these regions from observations. |Man-| 
jdelj ( |2010[ ) described the statistical procedure for inference 
on population parameters based on a limited set of uncer¬ 
tain observations. Farr et al. (20151 introduced a classifica¬ 
tion process based on models for subpopulation distributions 
whereby model parameters are fit as part of the classifica¬ 
tion process. Meanwhile, a number of unmodeled or weakly 
modeled schemes for extracting a distribution and classify¬ 
ing subpopulations exist, ranging from Dirichlet processes to 
k-means clustering (using the weighted posterior probability 
density functions of individual observations). The preferred 
algorithm depends on the priors that one wishes to place on 
the population distribution: modeled approaches require a 
choice of the shape (model) of the underlying distribution, 
which can aid measurement accuracy over unmodeled clas¬ 
sification schemes at the expense of introducing bias if the 
model shape does not match the observed population. In 
general, however, clustering should be successful whenever 
the presence of a break in the distribution can be estab¬ 
lished. As discussed above, for the BH-NS - BH-BH division 
this requires several tens of observations of each binary type 
if the distribution is similar to that described in Sec. HlBy 
contrast, the NS-NS and BH-NS subpopulation are cleanly 
separated in chirp mass, i.e., the measurement uncertainty 
a is very small relative to the size of the break in the chirp- 
mass distribution. In this case, the existence of a break can 

be established as long as ^N- 


is at least a few, i.e.. 


10 observations are sufficient. 

Counting sourees. Even if the subpopulations are clas¬ 
sified and their boundaries are accurately determined, indi¬ 
vidual sources cannot necessarily be perfectly classified as 
belonging to a given subpopulation. Nonetheless, the num¬ 
ber of sources in each subpopulation can be measured with 
some accuracy (Farr et al. 20151. The total mis-classified 
source fraction due to the uncertainty in cluster boundaries 

is < + -^A”^ — which, for the typical parameters 

we have considered here, is less than the fractional Poisson 
fluctuation in the subpopulation counts, ~ Therefore, 
we expect estimates of merger rates for various binary types 
to be limited by Poisson counting statistics on the number 
of events rather than classihcation errors. 


5 SUMMARY AND FUTURE DIRECTIONS 

We considered a plausible distribution of masses of merging 
NS-NS, BH-NS, and BH-BH binaries, whose gravitational- 
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wave signatures would be detectable by advanced ground- 
based gravitational-wave interferometers LIGO and Virgo. 
We found that NS-NS and BH-NS binaries were clearly sep¬ 
arated in chirp mass, while BH-NS and BH-BH distributions 
had some overlap in chirp mass, but were widely separated in 
the mass ratio. We evaluated the measurement accuracy of 
mass parameters using the LALInf erence toolkit that will be 
used for parameter estimation on LIGO-Virgo candidates, 
considering a variety of spin configurations. We concluded 
that if the mass distributions are known in advance, there 
will be no confusion on the NS-NS - BH-NS interface, and 
less than 1% of binaries could be incorrectly classihed on 
the BH-NS - BH-BH interface. We then considered the case 
where mass distributions are not known in advance, and 
determined that a few tens of detections are sufficient to 
correctly cluster the subpopulations in the two-dimensional 
mass space if the underlying population is similar to the 
standard model of Dominik et al. (20141. Although individ¬ 
ual sources may still be mis-classified, we anticipate that 
this will not be a signihcant source of error in the estimates 
of merger rates for various subpopulations. 

We considered a population synthesis model in which 
the supernova prescription, particularly the amount of fall¬ 
back on a newly formed BH, is designed to reproduce the ap¬ 
parent observational mass gap between NS and BH masses. 
If the NS and BH mass distributions are instead assumed 
to be continuous, as in the delayed supernova model of|Bel-| 


czynski et al. (20121, the quantitative results change: NS- 


NS and BH-NS chirp mass distributions overlap (Dominik 
et al.|20f4 1. However, the first percentile of the BH-NS chirp 


mass distribution is still significantly larger than the 99th 
percentile of the NS-NS chirp mass distribution, and simi¬ 
larly for the second and 98th percentiles of the BH-BH and 
BH-NS chirp mass distributions, respectively. In fact, the 
subpopulations are sufficiently concentrated away from the 
boundaries that the conclusions we have reached remain ro¬ 
bust, and clustering for the delayed supernova model is still 
possible with tens of observations. 

Another complication could be the low number of ob¬ 
servations of a particular species of binaries; for example, in 
the high BH supernova kick model of |Dominik et al!] ( |2014[ ), 
rates of BH-NS binaries are very strongly suppressed rela¬ 
tive to the other binary types. In this case, the number of 
BH-NS detections may be too small to identify a distinct 
cluster and appropriately classify them in the absence of a 
priori distribution models. 

This Letter shows that, at least for the physical mod¬ 
els considered, it should be possible to cluster and classify 
populations of NS-NS, BH-NS, and BH-BH merging bina¬ 
ries based on their gravitational-wave signatures. However, 
further rapid progress is required to address a number of 
practical challenges given the short timescale to anticipated 


GW detections ( 

Abadie et al. 2010 LIGO Scientific Col- 

laboration et al. 

20131. Actual clustering and classification 


schemes need to be implemented and tested for robustness 
on a variety of plausible source populations, evaluating the 
relative trade-offs of modeled approaches which require a 
choice of the shape (model) of the underlying distribution 
(Farr et al. 20151 vs. unmodeled classihcation schemes which 
are less sensitive to a priori assumptions about the distri¬ 
bution shape. Subpopulation rate estimates must take se¬ 
lection biases into account (Berry et al.|2014 Mandel et al. 


20151. The impact of other sources of parameter-estimation 


errors (e.g., systematics due to imperfect waveform knowl¬ 
edge, Aasi et ^|2013 Favata|[2014 1 needs to be evaluated. 
Classification should also be extended to consider other pa¬ 
rameters beyond component masses, such as spins and spin- 
orbit misalignment angles (e.g.. Vitale et al.||2014 l, as well 
as additional subpopulations, such as dynamically formed 
binaries in dense stellar environments. 
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